Exact solution of the Sutherland model 
with arbitrary internal symmetry 



Yusuke Kato and Yoshio Kuramoto 
Department of Physics, Tohoku University, 
Sendai 980, Japan 

Abstract 

An elementary theory is presented for solving the Sutherland model with ar- 
bitrary internal symmetry such as SU(z/) or a supersymmetry SU(V, [i). The 
ground-state wave function and all the energy levels are derived. One starts 
with solving a variant of the model with distinguishable particles, and then 
(anti)symmetrizes the solution. The theory is also applied to various lattice 
versions of the model. It is proved that the Gutzwiller-type wave function 
is not only an eigenstate of the supersymmetric t-J model, but is indeed the 
ground state. 
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There has been an increasing interest in a family of one-dimensional models with inter- 
action proportional to inverse square of the distance. With the periodic boundary condition, 
the model is called the Sutherland model Several years ago, Haldane |2| and Shastry || 
found that an S — 1/2 spin chain with long-ranged exchange (HS model) is solvable and 
that it is a lattice version of the original Sutherland model. After their works, the charge 
degree of freedom, hole, was introduced into the HS model and long-ranged supersymmetric 
t- J model has been proposed by Kuramoto and Yokoyama [§J. Ever after, the Sutherland 
model, HS model, and the t-J model were generalized to multi-component systems and the 
generalized models have been investigated intensively |5|-|7| • in view of fundamental physical 
interest in exact spin and charge dynamics in the t-Jmodel, for instance, it is highly desirable 
to develop a simple and flexible theoretical scheme for systems with internal symmetries. 

Sutherland and Shastry applied the asymptotic Bethe ansatz (ABA) to the multi- 
component Sutherland model and derived the spectrum and degeneracy, and discussed the 
thermodynamics ||. In spite of its power and simplicity, however, there are still unsolved 
issues in the ABA. One is about the validity of the asymptotic region in the system with 
a finite size. Another is about the ill-defined phase shift for equal momenta with different 
spins. Hence alternative approach is well worth a trial. 

In this paper, we present an elementary theory for deriving all the energy levels of the 
Sutherland model with arbitrary internal symmetry. The main idea of our method is firstly 
to generalize the calculation in ref. to a modified model for distinguishable particles. 
After that, we take into account the internal degree of freedom of particles and symmetrize 
or antisymmetrize the wave function to represent identical particles. The present theory 
has the following outstanding features. Firstly, our method identifies easily the ground 
state and gives the explicit form of the wave function for the most general models with 
the SU^v, fi) supersymmetry. Secondly, it gives a microscopic derivation of the energy as a 
functional of the momentum distribution function. Such an expression with correct account 
of the degeneracy is important especially for investigation of thermodynamics. Finally, it is 
also applicable to lattice models such as the multi-component t-J model. Because of these 
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features, this theory should provide a basis for further development such as investigation of 
dynamical properties for systems with internal symmetries. 
We consider the following model: 



a 2 2tt 2 A(A-CMg) 



(1) 



where is the exchange operator of coordinates of particles i and j ||. The size of the 
system is given by L, £ is 1 or —1 and the dimensionless coupling parameter A is positive. 
Since TC does not depend on the internal symmetry of particles, we consider only the orbital 
part of the wave function ^ ({x}). For the moment, we regard particles as distinguishable 
and does not impose any permutation properties on the wave functions. Here we write \1/ 
in the form \I/ = ^o.c^ where ^o,c is the "absolute ground state", by which we mean the 
ground state in the case where no restrictions to the symmetry of the eigenfunction are 
imposed. For \J/ ,c we can replace (M^ in Eq.(|l]) by 1 since the wave function satisfying 
(Mij^o^ = \I/o,c minimizes the repulsion. Then the result of ref. [l[ for the single component 
system can be used, and explicit form of \l/o,c i s given by 
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with the eigenenergy 
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Here N is the total number of particles. 

Note that \I/o,c i s the totally symmetric (antisymmetric) function when ( = 1 (—1). 
Considering $ as a function of the variables Zj = exp (i2irXj/L), we obtain the eigenvalue 
equation for $: 
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The eigenvalue e is related with the original eigenenergy E by e = [L/(2tt)} 2 (E — Eq) . 

Since the system has the translational symmetry, we consider eigenfunctions with the 
total momentum Q. We take a complete basis set of plane waves 

N 

i i 

where P — {p (1) ,p (2) , • • • ,p (N)) is an N-th order permutation and k = («i, k 2 , ■ ■ ■ , kat) 
is a set of integers ordered so that Ki > k 2 > • • • > k n . The action of W to Kj p gives 

H'<t> K ,p = £ [k? + A (iV - 1) /ti k ,p + 2A ^ h pU)p{k) 4>^p. (7) 

After some calculation, we obtain 

h p (j) P (k)4>K,p = -Kk4>K,p + R (\ z p(t)}) > ( 8 ) 

where i? is given by 

^|y-*k-i ( K . _ Kk _ ^ (zp^z"^) KjP (if > K fe + 2) , 



i2 



(9) 



(otherwise) . 

We then define the order of the basis. Let k' = (k'i, k' 2 , ■ ■ ■ , k' n ) be another set of 
momenta. We write k' < k if the first nonvanishing difference positive. For 

instance, (1, 1, 1, 1, 1) < (2, 1, 1, 1, 0) < (2, 2, 1, 0, 0). Similarly we define the order of the 
permutation P = (p(l) , p(2) , ■ ■ ■ ,p(n)). If P 1 is another permutation, we write P' < P if the 
first nonvanishing difference p{i)—p'{i) is positive. Each function of the basis is characterized 
by k and the permutation P. We write (/«', P') < (k, P) if k! < k, or k' = k and P' < P. 
The off-diagonal elements of T~C' come from the second term of Eq.([7p. From Eq.(|9|) we see 
that (V, P'lTi'lK, P) is zero when (k', P') > (k, P). Writing out 7i' in the ordered basis, the 
matrix is an upper triangular one in which all matrix elements below the diagonal vanish. 
We can obtain all the eigenenergies from the diagonal elements which come from the first 
terms of the right hand side of Eqs.(^) and (§) as follows: 

e K , P = e K = E ( K ? + A (JV - 1) Kj ) - 2A £ «* 

3 j<k 

= Eh 2 + A(iV + l-2j)s]. (10) 



Here we introduce the momentum distribution function v{k) by v (k) = J2j $ k) with 
S k) the Kronecker's delta symbol. Then we can rewrite the expression (|T0| ) as a func- 
tional of v (k): 

oo \ 

e K = £ tfi/(A0 + ££ !*-*>(*)*'(*')■ ( U ) 

fc=— oo k,k' 

All the eigenf unctions are written in the form: 

^K.P = 4>K,P + J2Yl a K',P'4>K',P' • (12) 
K'<K P> 

Each eigenvalue depends only on k and is independent of permutation P. Thus, there is 
one-to-one correspondence between each eigenstate $ Kj p of the present model and that of 
the free system (A = 0) given by (f) Kt p. 

Now we consider particles as SU(v) fermions and calculate the ground state under the 
given color distribution {M a }^ =1 where M a is the total number of particles with the color a. 
Firstly we consider the case £ = 1. Since the absolute ground state \l/o,i is totally symmetric, 
the momenta of particles with the same color are ordered non-equally (• • • Kj > Kj + i ■ ■ •). We 
introduce the momentum distribution of particles with the color o by 

(k) = ^25(K j ,k)5(a,cr j ). (13) 
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In this case v a (k) is either or 1. The energy relative to E is given by Eq.(|ll|) with 
v (k) = J^a^a (k). Let us consider the following momentum distribution: 

u a (k) = e(M a /2-\k\), (14) 

where 9 (k) is the step function, and each M a is taken to be odd. It is clear that Eq.(|T4"D 
gives the minimum for both the first and second terms in Eq.flTTD as SU{y) fermions. Hence 
this distribution gives the ground state. 

Let Kg be the set of the ordered momenta corresponding to the momentum dis- 
tribution (0). The ground-state wave function <3> g is given by antisymmetrization of 
$K g ,pP ({Mr}) X ({°j}) where x is a color function and P ({(Jj}) is given by 
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P(M) = US\M ^5(a jl a)\ . (15) 

Firstly we antisymmetrize the coordinate of particles with the same color. Because of the 
compact momentum distribution, the only part in $k 9 ,p that survives antisymmetrization is 
<pK g ,p- Using the formula for the Vandermonde determinant, the antisymmetrization gives 

e (-i) p n ^ +1)/2+i = n ^ Mv ~ m n & - *o . m 

P j=l j=l j<k 

where (— 1) P denotes the sign of the permutation. We must multiply this by the antisym- 
metric color function for different colors. The result is given by 

jv 
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2-sgn {Gj - a k ) 



(17) 



Ha and Haldane proved that this wave function is an eigenfunction of a model equivalent 
to Eq.(^) with £ = 1, and conjectured that it provides the ground-state wave function ||. 
We have proved here that their conjecture is correct. To see the equivalence between the 
models we note that —My in Eq.(^) can be replaced by the color exchange operator for 
SU (v) fermions. On the other hand if M a is even the ground state is degenerate because 
the highest momentum occupied can be either with plus or minus sign. Provided that N/v 
is an odd integer, the minimum energy distribution among all the color distribution {M a } 
is given by M a = N/ v for all cr's. 

Next we consider the case ( = —1 for SU (is) fermions. In this case the absolute ground 
state ^0,-1 is totally antisymmetric and the momentum distribution is bosonic: v a (k) = 
0,1,2 The ground state is given by v a (k) = M a 5(k,0) for all cr's and the ground- 
state energy is just E . The corresponding eigenfunction is ^ = ^o,-i ({%}) P (Wj})- Any 
distribution {M a } leads to the same energy given by Eq.([|). Hence there is a degeneracy 
(N + u- 1)1/ [{v - 1)\N\] in the ground state. 

In the cases where particles are SU(v) bosons or the mixture of bosons and fermions, the 
model given by Eq.([TJ) can be solved in a similar way. Thus we consider the most general 
SU(is,fi) model, i.e., particles consisting of SU(p) fermions and SU(is) bosons. The wave 
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function for the ground state has not yet been reported for this model. In this case, we can 
rewrite as 

{— Pi a (if both i- and 7-th particles are fermions), 
(18) 
Pij (otherwise) . 

We first consider the case £ = 1. The absolute ground state is given by 

N B N F N B N F 

*o,i = n \tj - &i A n - n n & - ( i9 ) 

J<fe j<fe J k 

where £j (a>j) and <7 B are variables for the complex coordinate and the color of j-th 
boson (fermion), respectively, and N B (N F ) is the total number of bosons (fermions). 

Let the momentum distribution for bosons (fermions) with the color a ((3) be given by 
{"a ( k )} ({ Z/ J( fc )})> and the color distribution by {M^}^ NMj}^_J where each 
Ml is taken to be odd. Since the absolute ground state fll9l) is symmetric, z/ B (fc)'s take all 
non-negative integers : 0, 1, 2, ■ ■ • and v% (fc)'s do or 1. The minimum energy distribution 
is given by 

v*(k) = M*8(k,0), 4{k) = 6(M*/2-\k\). (20) 
The corresponding wave function \1/ is given by ^o.i'&g where 



x exp ^sgn (aj - P ({m b }) P ({mJ}) , (21) 

with K(a) = (Mj' 6 — l)/2. Similarly, the ground state in the case ( = —1 is given by 
^o,-i^p where 

*;=ff«; i ' ( " ?) ffte-&) i( "' v?) 

x exp [^sgn (af - af)] P ({M Q B }) P ({m|}) . (22) 

Our method is also applicable to lattice models by taking the limit A — > oo |||1(J. In the 
strong coupling limit, a part of the potential A 2 J2j<k sin -2 [tc (xj — Xk) / L] enforces particles 
to localize with a lattice spacing L/N. Up to order O (A), the Hamiltonian decouples into 



T~t = 7~Lci& + 'Hut + E 



Mad ■ 



(23) 



where TC e ia,, Hint, an d -^Mad are the elastic and lattice Hamiltonians and the Madelung energy, 
respectively. We define Uj as the displacement from the j-th lattice point (j = 1, 2, • • • , N) 
and introduce a function 

J (j ) = 3sin- 4 (^)-2sin- 2 (^). (24) 

Then we obtain 



E d(j-k) 
.*(#) 



II,: 



4A2 (t) EdU-k)u jUk , (25) 



L 



j<k 



where j and k are no longer indices of particles but those of sites. We rewrite H, e i& as 

N-l /„n 2 



(26) 



where ¥ g [b q J is the creation (annihilation) operator of the phonon field. In deriving the 
expression (|26"D, we have used the following result: 

iV 4 -l 
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where Pjk is the color exchange operator between the j-th and k-th sites. When ((, u, fi) = 
(1, 0, 2) or (—1, 2, 0), Ti.\ a t becomes the Hamiltonian of the HS model. In the case of ((, z/, fi) = 
(1, 1, 2) or (—1, 2, 1), it becomes the Hamitonian of the supersymmetric t-J model. 

We now derive the wave function = ^l^g of the ground state with ( = 1 in the 
limit of A — > oo. The absolute ground state \l/ ,i becomes 

2 JV-1 



#0,1 ex P 



Mt) T, ( l(N-q)u(q)u(N-q) 



q=0 



x (numerical factor) 



(30) 
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where u (q) is the Fourier transform of the displacement Uj. The \l/o,i corresponds to the zero- 
point motion of the phonon field. Factoring out the phonon part, we obtain the ground state 
of the SUiy^jj) lattice model, which is nothing but $ 9 given by Eq.(^TJ). The remaining 
effect of \l/o,i is the exclusion of multiple occupation of each site. In the special case of 
SU(1, 2), it is therefore proved that the Gutzwiller wave function is indeed the ground state 
of the t-J model [|. 

Thus the present theory gives a unified treatment of the family of Sutherland-type mod- 
els. In every case, our theory derives all the energy levels by a simple calculation. We have 
proved in this way the previous conjectures on the ground states. Furthermore simplicity of 
the method permits us to derive new results for the ground-state wave function of SU (v, /x) 
models in both the continuum space and the lattice. 

Lastly, we make a brief remark on another application of our theory. Similar approach 
is useful to the Calogero model with the internal symmetry as well []12" |. As a matter 
of fact, the Calogero model is much more tractable than the Sutherland model. For the 
Calogero model, the local operators for energy boost were already found P,|TT||. By using 
such operators, we have obtained all the eigenstates of the SU (u) Calogero model and proved 
that the eigenfunction proposed in refs. ]7|jr2(] is the exact ground state. The details will be 
discussed elsewhere. 
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